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Abstract 

We study the energy released from phase-transition induced collapse of neu- 
^ ! tron stars, which results in large amplitude stellar oscillations. To model this 

^vq I process we use a Newtonian hydrodynamic code, with a high resolution shock- 

ed ■ capturing scheme. The physical process considered is a sudden phase transition 

from normal nuclear matter to a mixed phase of quark and nuclear matter. We 
*0 I show that both the temperature and the density at the neutrinosphere oscil- 

^-^ • late with time. However, they are nearly 180° out of phase. Consequently, 

extremely intense, pulsating neutrino/antineutrino and leptonic pair fluxes will 
be emitted. During this stage several mass ejecta can be ejected from the stel- 
lar surface by the neutrinos and antineutrinos. These ejecta can be further 
accelerated to relativistic speeds by the electron/positron pairs, created by the 
5^ I neutrino and antineutrino annihilation outside the stellar surface. We suggest 

that this process may be a possible mechanism for short Gamma-Ray Bursts. 

1 Introduction 

Recently, by using simulations performed with the Newtonian numerical code in- 
troduced in P, it was shown that the resulting quark star, produced by the phase 
transition induced collapse of a neutron star, will undergo a series of oscillations. The 
collapse process with a conformally fiat approximation to general relativity was also 



>$ 



simulated in [2] . The works of [H [2] focus on the gravitational wave signals emitted 
by the collapse process. 

It is the purpose of the present paper to consider another important implication 
of this result, namely, the effect of the oscillations of the newly formed quark star 
on the neutrino emission. The oscillations can enhance the neutrino emission rate 
in a pulsating manner, and the neutrinos are emitted in a much shorter time scale. 
Moreover, through the process of neutrino-antineutrino annihilation, a large amount 
of electron-positron pairs is also produced. Therefore the neutron-quark phase tran- 
sition in compact objects may be the energy source of GRBs [3]. Such a model 
can also explain the lack of detection of a neutron star or pulsar formed in the SN 
1987A |1], by assuming that the newly formed neutron star at the center of SN 1987A 
underwent a phase transition after the neutrino trapping timescale (~ 10 s). Conse- 
quently, the compact remnant of SN 1987A may be a strange quark star, which has a 
softer equation of state than that of neutron star matter. Such a phase transition can 
induce stellar collapse and result in large amplitude stellar oscillations. Extremely in- 
tense pulsating neutrino fluxes, with submillisecond period and with neutrino energy 
(greater than 30 MeV), can be emitted because the oscillations of the temperature 
and density are out of phase almost 180 °. If this is true, the current X-ray emis- 
sion from the compact remnant of SN 1987A will be lower than 10^^ erg s~^, and it 
should be a thermal bremsstrahlung spectrum for a bare strange star with a surface 
temperature of around ~ 107 K. 

2 Description of the phase transition and of the 
equations of state 

The numerical code is based on the three-dimensional numerical simulation in New- 
tonian hydrodynamics and gravity. The quark matter of the mixed phase is described 
by the MIT bag model and the normal nuclear matter is described by an ideal fluid 
EOS. This code has been used to study the gravitational wave emission from the 
phase-induced collapse of the neutron stars P, where a detailed discussion can be 
found. In the simulations, we introduced a very low density atmosphere outside the 
star. The "artiflcial" atmosphere is not physical, but it is important for the stability 
of the hydrodynamical code. This is due to the problem that the hydrodynamical 
codes cannot in general handle vacuum regions where the density is zero. In order to 
avoid a signiflcant influence of the atmosphere on the dynamics of the physical sys- 
tem, it is necessary to choose the density of the atmosphere patm to be much smaller 
than the density scale of interest. For the results reported in this paper, we set patm 
to be 3 X 10^ g/cm . 

In our simulations, we will not simulate the phase transition process. Instead, 
we assume that a fast phase transition has happened (e.g., via a detonation mode) 



so that the initial neutron star has converted to a quark star in a timescale shorter 
than the dynamical timescale of the system. We assume that the normal matter 
inside the initial neutron star has suddenly changed to quark matter at t = 0. This 
is achieved by changing the EOS at t = after the initial hydrostatic equilibrium 
neutron star has been constructed. We then simulate the resulting dynamics of the 
system triggered by the collapse. 

The equation of state (EOS) for neutron stars is highly uncertain. We could 
try all possible existing realistic EOS in our study. However, the main purpose of 
this paper is to demonstrate that during the phase-transition induced collapse of 
a neutron star extremely intense, pulsating and very high energy neutrinos can be 
emitted. The effect is governed mainly by the amount of pressure reduction after the 
phase transition as compared to the initial neutron star model. For simplicity we will 
use a polytropic EOS for the initial equilibrium neutron star: 

P = kop'^^ (1) 

where fcg and Fq are constants. On the initial time slice, we also need to specify the 
specific internal energy e. For the polytropic EOS, the thermodynamically consistent 
e is given by 

pr«-\ (2) 



Fo-1 
Note that the pressure in Eq. ([T]) can also be written as 

P = (Fo - l)pe. (3) 

In this paper we use a mixed phase EOS to mimic a quark star core covered by 
normal nuclear matter. This EOS model consists of two parts: (i) a mixed phase of 
quark and nuclear matter in the core at density higher than a certain critical value 
Ptr (quark seeds can spontaneously produce everywhere when p > ptr) (ii) a normal 
nuclear matter region extending from p < ptr to the surface of the star. A more 
detailed discussion about such hybrid quark stars can be found in [T]. Explicitly, the 
pressure is given by 

{aPq + (1 - a)Pn, for p > ptr 
(4) 
Pn, for p < Ptr, 

where 

P^ = \{p + pe-AB), (5) 

is the pressure contribution of the quark matter, 

Pn = (F„ - l)pe, (6) 



and 

f (p-Ptr)/(Pg-Ptr), for Ptr < p < pq, 

(7) 

[ 1, for Pq < P, 

is defined to be the scale factor of tlie mixed phase p-j. r„ is not necessarily equal to 
Fq. As for the quark matter we have assumed that it is described by the MIT bag 
model. It should be noticed that Pg is not in the usual form of P = (ptot — 4-B)/3, 
where ptot is the (rest frame) total energy density, and B is the bag constant. It is 
because in our Newtonian simulations, we use the rest mass density p and specific 
internal energy e as fundamental variables in the hydrodynamics equations. 

The total energy density ptot, which includes the rest mass contribution, is decom- 
posed as Ptot = P + P^- We choose r„ < Fq in our simulations to take into account 
the possibility that the nuclear matter may not be stable during the phase transition 
process, and hence some quark seeds could appear inside the nuclear matter, or the 
convection, which can occur during the phase transition process, can mix some quark 
matter with the nuclear matter. In the presence of the quark seeds in the nuclear 
matter, the effective adiabatic index will be reduced. The possible values of B^^'^ 
range from 145 MeV to 190 MeV El El [H] • For p > pg, the quarks will be decon- 
fined from nucleons. The value of Pg is model dependent; it could range from 4 to 8 
Pnuc [El m [10], where pnuc = 2.8 x 10^^ g cm~^ is the nuclear density. 

The properties of a mixed quark-hadron phase and its implications for hybrid star 
structure were considered in [TT]. In this study it was assumed that charge neutrality 
must be enforced only as an overall constraint, an not separately on each of the phases. 
In the mixed phase the fraction x of the volume occupied by the quark phase can 
be taken as the basic independent variable, being more convenient than the baryon 
density. The pressures of the two phases in equilibrium are equal, but vary with the 
quark fraction. Once the properties of the two phases are known, with the use of the 
relationships between the quark chemical potentials and the independent conserved 
quantities like baryon number and electron charge, the equation of state of the matter 
in the mixed phase can be obtained. In Fig. [1] we present the equation of state of 
the mixed phase obtained in [11] for a compression modulus K = 240 MeV and an 
effective mass at saturation density of m*/m = 0.78, and the equation of state Eq. (jl]) 
used in the present simulation, respectively. 

In the simulations, we choose Fq = 2, F„ = 1.85, B^^^ = 160 MeV and p, = 9pnuc- 
The transition density ptr is defined to be at the point where Pq vanishes initially. 




log(p) 



Figure 1: Comparison of the equation of state of the mixed quark-hadron phase 
proposed in [TT] for K = 240 MeV and effective mass m*/m = 0.78 (dotted curve), 
and Eq. (jlj), the equation of state if the mixed phase used in the present simulations 
(sohd curve). 



3 Emission of neutrinos and e^ pairs 



The neutrino luminosity is given by 



4:7ir^c ■ 



Eu d^Pp 



27ih J I + exp{Ei, — jj,i,/kT^ 



(8) 



where fi,y is the neutrino chemical potential. Taking fi^ = 0, the neutrino luminosity 
emitted from the neutrinosphere is given by L^, = TnRlacT^ /16, where a = 4a /c is the 
radiation constant, a is the Stefan-Boltzmann constant, and T,y is the temperature of 
the neutrinosphere. If we assume equal luminosities for neutrinos and antineutrinos, 
the combined luminosity for a single neutrino flavor is 



Lj,.77 = L^ + Ltj= -7rR,McT„ 



(9) 



The effect of coherent forward scattering must be taken into account when consid- 
ering the oscillations of neutrinos traveling through matter [15]. Although different 
flavor neutrinos have different R,^, yet they have approximately the same value of lu- 
minosity for all flavors [121 [IT]. Therefore the total luminosity is around three times 
of a single neutrino flavor luminosity 



Lu^.v. + -^i/„,77„ + -^z/..77. = --iTRlacT, 



(10) 



Using Ri, and T^, obtained from the last Section, we compute the the neutrino 
luminosity as a function of time. The results are shown in Fig. |2j 
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Figure 2: Neutrino luminosity versus time for M = 1.55Mq, (left upper figure), 
M = 1.7Mq, (right upper figure), M = l.SM©, (left lower figure), M = I.QMq, 
(right lower figure). 



The peak luminosities range from 10^^ to 10^^ ergs/s; the pulsating period of the 
luminosity is the same as that of the temperature and of the density. 

Neutrinos and antineutrinos can become electron and positron pairs via the neu- 
trino and antineutrino annihilation process (i/ + z7 — )■ e~ + e"*"). The total neutrino 
and antineutrino annihilation rate can be given as follows [131 lE^ 



^uv- 



'-^j^D[kUt)r £eir)4nr^dr 



7GlDTT^Ci5) 8tt'^ 
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Rl{kT,)\ 



(11) 
(12) 



where 9(r) = 27r^(l — xY{x^ + 4x + 5)/3, x = J\ — Rl/r"^, Ty(t) is the temperature 
at the neutrinosphere at time t, G\ = 5.29 x 10"^'* is the Fermi constant, ( is the 
Riemann zeta function, and D is a numerical value depending on the pair creation 
processes (e.g. experimental results indicate that Di = 1.23 for z/g Ue and D2 = 0.814 
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for z/^ Vjx and z/^ z/-). To obtain the total neutrino annihilation rate from all species, 
^e t'e; ^fj. ^ji ^nd z/^ z/-, we sum up the energy rate for each single flavor, 



28G2,7r6C(5) 



9c5/i6 



{D^ + 2D,)Rt{knr. 



(13) 



Fig. [3] shows the rate of energy carried away by the electron/positron pairs pro- 
duced through neutrino annihilation, which varies from ~ lO^^ergs/s to ~ lO^^ergs/s. 
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Figure 3: Electron-positron luminosity versus time for M = 1.55Mq, (left upper 
figure), M = 1.7Mq, (right upper figure), M = l.SM©, (left lower figure), M = 
1.9Mq, (right lower figure). 



It is interesting to note that almost all neutrinos can be annihilated into electron- 
positron pairs at the peak because of the extremely high density and high energy of 
the neutrinos. In particular the rest mass of the electrons/positrons is much smaller 
than kTy. 



4 Discussions and final remarks 

In this paper we have studied the possible consequences of the phase-induced collapse 
of neutron stars to strange stars. We have found that both the density and the 
temperature inside the star will oscillate with the same period, but almost 180° out 
of phase, which will result in the emission of intense pulsating neutrinos and pairs. 
We want to point out that the intense pulse neutrino luminosity can be maintained 
due to the oscillatory fluid motion, which can carry thermal energy directly from the 
stellar core to the surface. This process can replenish the energy loss of neutrino 
emission much quicker than the neutrino diffusion process. A large fraction of the 
neutrino energy, roughly (1-1/e), will be absorbed by the matter very near the stellar 
surface. When this amount of energy exceeds the gravitational binding energy, some 
mass near the stellar surface will be ejected, and this mass will be further accelerated 
by absorbing pairs created from the neutrino and antineutrino annihilation processes 
outside the star. Although matter will be ejected periodically, each ejecta can have 
different masses and Lorentz factors, and therefore the intrinsic period could not be 
observed. We suggest that the collisions among these ejecta may produce short GRBs. 
Our numerical simulations are simulating a spherically symmetric and non-rotating 
collapsing stellar object, and they also do not contain magnetic fields. Therefore the 
radiation emission produced in this model is isotropic. However, a realistic neutron 
star should have finite angular momentum and strong magnetic field, and hence these 
two factors could produce asymmetric mass ejection. This effect will be considered 
in future work. 

The phase-transition from a neutron star to a strange star was simulated in [18j, 
with the conclusion that this process is most likely not a gamma-ray burst mechanism. 
They mimic the phase-transition by the arbitrary motion of a piston deep within the 
star, and they have found that the mechanic wave will eject ~ IO^^Mq baryons, which 
causes the baryon contamination for the gamma-ray bursts. In our simulations, we 
assume a sudden change of equation of state to mimic the phase-transition, and we 
use the Newtonian hydrodynamic code to study the response of the stellar interior 
after such a sudden change of the EOS. In our simulations we find that the mass 
ejection by the motion of the fluid is very small. We estimate that the major mass 
ejection would result from the heating of neutrinos and pairs on the stellar crust, 
which is not modeled in the simulations. Our total energy output and total mass in 
ejecta are close to that of [TB]. However, the neutrino energy injection is pulsating, 
and hence the mass ejection is also pulsating. The mass of individual ejecta range 
from ~ 10~^ to ~ IO^^Mq, with output energy in the range of 10^*^ to 10^'' ergs. 
Therefore, some ejecta cannot be relativistic, and they cannot contribute to GRBs. 
However, there are still many relativistic ejecta in each simulation model, which can 
have Lorentz factors >100, and with a total energy of ~ 10^° — 10^^ ergs. This could 
be a possible mechanism for short GRBs. 
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